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SUMMARY 

This Investigation deals with the stability characteristics of oil 
journal bearings. Including the effect of elastic distortions In the bearing 
liner. Graphical results are presented for (1) steady-state load, (2) 
stiffness and damping coefficients, and (3) the stability. These results are 
given for various slenderness ratios, eccentricity ratios, and elasticity 
parameters. The lubricant Is first assumed to be Isovlscous. The analysis Is 
then extended to the case of a pressure-dependent viscosity. It has been 
found that stability decreases with Increase of the elasticity parameter of 
the bearing liner for heavily loaded bearings. 

NOMENCLATURE 

a,b Inside and outside radius of the bearing 

B HqOwR^/C^ 

C 2 + \/jj 

c radial clearance 

D journal diameter 


♦National Research Council - NASA Research Associate. 



Drr, Dr«p, 
Dq>r , Dcfxp 


dm,n 

E 


e 

F 

G 


H 

h, h 
1 

Krr, Kr<p, 
K<pr , Kcptp 


L 

M, M 
m, n 
P, P 


Po* Pi ♦ P2 


q, q 

R 

x, y, z 
e, y, z 
T 
t 

u, v, w 
W, W 

a 

13 


dimensionless damping coefficients (0 = Dc 3 /n 0 R 3 L), where the 
first and second subscripts refer to the component of the force 
and the direction of the velocity respectively 

distortion of the m,n harmonic 

Young's modulus 

eccentricity 

elasticity parameter, F = n 0 R 3 w/c 3 E 
aE 

thickness of bearing liner 
film thickness, h = h/c 

yrr 

dimensionless stiffness coefficients (K = Kc 3 /n 0 uR 3 L) , where the 
first and second subscripts refer to the component of the force 
and the direction of the displacement respectively 

length of bearing 

mass of rotor, M, mass parameter (M = Mcu 2 /W) 
axial and circumferential harmonic 
pressure, (p = pc 2 /n 0 wR 2 ) 

steady-state and perturbed (dynamic) dimensionless 
pressures 

modified pressure, (q = 1/B(1 - e-P B )) 
journal radius 

circumferential, radial, and axial coordinates 
dimensionless coordinates, e = x/R, y = y/a, z = z/(L/2) 
nondlmenslonal time, T = w p t 
time 

radial, circumferential, and axial displacements 
steady-state load, (W = Wc 2 /no<*>R 3 L) 
piezo-viscosity coefficient 
harmonic phase angle 
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deformation of bearing surface 


6 

c, c q eccentricity ratio (e = e/c), e 0 steady-state eccentricity 

ratio 

n, n 0 absolute viscosity of oil, (n 0 = viscosity at Inlet condition) 

i|rj ,4>2 angular coordinates at which film commences and cavitates, 

respectively measured from minimum film thickness 

t stress 

©1 « + 

©2 ir + i ^2 

p Lame's constants 

<j Poisson's ratio 

<p, «p 0 attitude angle, <j> 0 = steady-state attitude angle 

a whirl ratio (Si = up/u) 

<a> angular velocity of journal 

wp angular velocity of whirl 

INTRODUCTION 

Theoretical research on flexible (soft shell) bearings with a rigid rotor 
was started with the work of Hlgglnson (Ref. 1) using a simplified method (the 
distortion Is proportional to the pressure). Since then many workers notably 
Hooke, Brighton, and O'Donoghue (Refs. 2 to 4), Conway and Lee (Refs. 5 
and 6), and Oh and Huebner (Ref. 7) solved the Journal bearing problem 
considering the effect of elastic distortions of the bearing liner. In 
Refs. (2) and (3) the analysis dealt with the Infinitely long bearing 
approximation. Brighton, Hooke, and O'Donoghue (Ref. 4) described the method 
of solution for finite bearings considering the effect of elastic 
distortions. This paper mentions some of the Important experimental findings 
of Carl (Ref. 8). Bozacl et al . (Ref. 9) gave theoretical and experimental 
Investigation of a finite journal bearing assuming that three-fourths of the 
bearing arc was a uniform thin ring and the remaining one-fourth was rigid. 
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Conway and Lee used long and short bearing approximations, respectively In 
their papers (Ref. 5 and 6). Oh and Huebner (Ref. 7) used a finite-element 
method which was applied to finite journal bearings. 

In recent years several workers have studied the effect of elastic 
distortion In dynamically loaded journal bearings (Refs. 10 to 13). From 
these studies the minimum film thickness as a function of crank angle of 
engine bearings can be determined. 

The purpose of the present study Is to Identify the threshold of oil 
whirl for a rigid rotor In self-acting oil journal bearings considering the 
elastic deformation of the bearing liner. The deformation of the liner Is 
estimated by solving the displacement components In the elastic medium. In 
the stability study both the translatory and rotating components of journal 
vibration have been considered. Using a first-order perturbation method the 
time-dependent Reynolds equation Is transformed In terms of steady-state and 
dynamic pressures. These pressures are solved by a finite-difference method 
using Reynolds boundary conditions. From the dynamic pressures the stiffness 
and damping coefficients are computed. These are then used In the equation of 
motion to find the mass parameter, which Is a measure of stability. In the 
absence of experimental data, the theoretical mass parameters obtained by this 
analysis have been compared with the available results of rigid bearings. 

THEORY 


Using the normal assumptions In the theory of hydrodynamic lubrication 
the differential equation In the bearing clearance of an oil-lubricated 
bearing as shown In Fig. 1 Is 


3 ,,3 l, 3 

M < h i? f M 


, k 3 . n/1 2 a®, ah ,ah. 

(h f z ) , 6 V R(1 - - at ) « * 12n„ <it> 


(') 
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Eq. (1) when nondlmenslonal 1 zed with the following substitutions. 



_z 

L 

2 


2 <>> 
p = PC 0 , T = w t , and Q = , will read as: 

n c P CO 

V R 


a_ 

ae 


(fi 3 ||) * (“> (S 3 % = 6(1 - 20 If) li * 120 If 

3 Z 3 Z 


( 2 ) 


The steady-state film thickness h In the case of a flexible bearing Is 

o 

given by 


fi o = 1 * c o cos 6 + c (e,z) ( 3 ) 

where 6 Is the deformation of the bearing surface. 

Before trying to find the solution of Eq. (1) satisfying the appropriate 
boundary conditions the elastic deformation 6 Is obtained In the following 
way: 

The method Is similar to that of Brighton et al. (Ref. 4). In the present 
calculation the three displacement components u, v, and w are solved 
simultaneously satisfying the boundary conditions. Brighton et al. (Ref. 4), 
however, used an approximate method to evaluate the displacements. The method 
Is briefly given below. 

The pressure distribution In the bearing clearance of the rigid bearing Is 
first calculated by solving the two-dimensional steady-state Reynolds equation. 
The film pressure Is then expressed In a double Fourier series of the form: 

I I 

P = E L P m ,„ '°s cos ' ne * < 4 > 

I 

where V" Indicates the first term of the series Is halved, and p 

/ > m,n 

and 0 are given by 
m, n 


5 


m, n wL 


2k J./2 


f ft 


p cos - cos ne dz do ( 


2ir L/2 


in 


p cos ^ - - sin ne dz de 


1/2 


(5) 


13 = tan 

m,n 


-i 


2ir J-/2 


[ i 


p cos - sin ne dz de 


2if X/2 


I I 


2m * z _ 
p cos — cos ne dz de 


The first term of the right-hand side of Eq. (4) Is 1/2 P q q . When 

the condition at the end of bearing (p = 0 at z = L/2) Is used, one can 

obtain p . Although this term does not contribute any deformation at 

0,0 

z = L/2, Its effect for other values of z Is already Included In the total 
deformation. The boundary conditions of the Inner radius are 

T rr = - p 

’ r6 - 0 I < 6 > 

T =0 
rz 

The outer surface Is rigidly enclosed by the housing preventing any 
displacement of the outer surface. The ends of the bearing are prevented from 
expanding axially, but are free to move circumferentially or radially. 

The displacement components In r, 6, and z directions were found 
from the pressure distribution which had been expressed In a Fourier series. 

These displacements were substituted In the stress-strain relationships 
using Lame's constants. The six components of stresses were then used In the 
equations of equilibrium to obtain the following three displacement equations. 
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2 2 
r d u C du tr n 2 N u /r , . n dv 

C — r + - (C + n ) + (C - 1 ) - 

dy y dy y y dy 


- (C ♦ l)^r v-^u + (C-l)^ — = 0 
y * a dy 


(1 + Cn 2 ) 

dy y dy y 


- (C - 1) 1 dM 
y dy 


(c + l)^ - ^ (c - i) r s 0 

y a y 


d w + 1 dw _ p 2 w_ 
dy 2 y dy y 2 


C^w - (C-l) - (C - 1) ^ - n 2 (C - 1) - = 0 

a dy y y 


( 7 ) 


J 


where 


2 + - , k 
v 


2 m ir a 
L 


The boundary conditions are: 


At y = 1, C 


^ . - i - (C - 2)(fl-S * a * L 2 w) 
dy y y a 2 


and 


dv u v 

— a - + - 


dy 

dw 

dy 


dw 

— = u 


( 8 ) 


at y=~,u = v= w = 0 

a 

Equation (7) was solved using finite-difference methods for unit 

pressure, y and R, and the values of d were obtained and expressed as 

m,n 


d = ^ 

m,n Rp 

The deformation A of the bearing surface will be 


(9) 


2m ir z 


A = V V P m n d m n cos ~ (cos ne + R ) 

m,n m,n L m,n 


( 10 ) 
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Knowing the distortion coefficients d m>n and using the expressions for 
p m>n and from Eq. (5) d at any point (e,z) was computed. 

Using p = pc 2 /n 0 «*)R 2 and z = z/(L/2), the radial deformation In the 
Inner bearing surface Is 


- 2(1 £ Vn Vn cos cos (n9 * 0 n..n )9(6 ’ 5> (11 > 


3 3 

Here p Is replaced by E/2 ( 1 + a) and F = n Q «R /Ec . 

Let us now solve Eq. (2). The boundary conditions are: 


and 


P(e,±l) = 0 

^(e,0) = 0 

az 

p(e,z) = 0 


»£ _ M _ o 
- ~ ae 
az 


( 12 ) 


- 6 - e i 


where e Is the angular coordinate where film cavitates. 

Assuming that the journal whirls about Its mean steady-state-position 

given by e and «p , for the first-order perturbation, the pressure 
oo 


and film thickness can be expressed as 

IT 


IT 


and 


where 


P = p o + p i c i e + p 2 c o*l e 


fi = fi + c.e^ T cos 6 + c 6. e^ sin 6 
0 1 0 1 


(13) 


e = e Q + c^e 


IT 


( 14 ) 


V = <P 0 + <P-| e 


IT 


Substituting Eq. (13) Into Eq. (2) and collecting the zeroth and the first- 

order terms for e, and e <p, , the following set of equations result 

1 o 1 
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( 15 ) 


• >b( 


+ 3(f) 2 ^ (fi o 2 ^ cos ej* - 
az \ az / 


cos eU -6 sin e + 1(12 q cos e) (16) 


and 


a_ 

ae 




♦ 3( 


t >2 ti (*° 2 f •)■ 


6 cos e + 1 ( 24Q sin e) - 1 


( 12a 
c 0 C 36 ) 


(17) 


Steady-state Solution 

The bearing pressure was first obtained from Eq. (15) assuming a constant 
film shape (l.e., h Q = 1 + c q cos 0) and using a finite-difference method 
(Gauss-Seldel ) with successive overtaxation scheme. The distribution was 
expressed as a double Fourier series as given by Eq. (4). The deformation 4/c 
(Eq. (11)) was then calculated for a given F using distortion coefficients 
from Eq. (9). The film thickness equation was then modified using Eq. (3). 

The process was repeated until a compatible film shape and pressure 
distribution was determined. The steady-state load and attitude angle are 
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Solution Under Dynamic Condition 

Having obtained the steady-state pressure distribution of a journal 
bearing having various L/D (0.5, 1.0, and 2.0), c q (0.6, 0.8, and 0.85), 
and F (0, 0.05, 0.1, 0.2, and 0.4) for a constant value of H/R = 0.3 and 
a = 0.4, the pressure distribution under dynamic condition was obtained. 
Liner thickness to radius ratio of 0.3 and Poisson's ratio of 0.4 were taken 
In the present calculation for the purpose of comparison of steady-state 
results of Brighton et al. (Ref. 4). The method, however, can handle any 
reasonable value of H/R and a up to 0.49. For this purpose Eqs. (16) 
and (17) are solved satisfying the appropriate boundary conditions by a 
finite-difference method using a successive overrelaxatlon scheme. 

Stiffness and Damping Coefficients 

Nondlmenslonal stiffness and damping coefficients can be shown to be 
(Ref. 15) 
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cos e de dz 



P 2 sin e de dz 


P 2 cos e de dz 




sin e de dz 


rr 



cos e de dz 


0 

w 



P 2 sin e de dz 


Q 





Pi 

o~ 


sin e de dz 


( 20 ) 
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Stability Characteristics 


Referring to Fig. 1, the equations of motion of the rigid journal, along 
the line of centers and Its perpendicular direction can be written as 


and 


M c 


a 2 H* 2 

d c_ , 0<i>' 


dt 


2 cl dt' 


= F r + W cos «p 


" c ■« H - 2 <if) <dt> 

dt 


= F - W sin «p 
<p 


( 21 ) 


( 22 ) 


F and F are the resultant film forces In the r and <p directions, 
r v 

For steady-state 


F + W cos <p = 0 
ro o 


(23) 


F - W sin <p = 0 \ 

«po o J 

Substituting Eqs. (14) and (23) Into Eqs. (21) and (22) and nondlmenslonallzlng, 
after neglecting second-order terms 

(-MQ 2 W + 1QD„ + K ) Cl + (1c QD r + e R + W sin <o ) <p, = 0 (24) 

rr rr 1 o r<p o r<p o 1 

(1QD + K ) e, + (-MQ 2 e W + IQ e D + c K + W cos » ) <p, = 0 (25) 

<pr <pr 1 ' o o qxp o <p<p o 1 

For a nontrivial solution the determinant of Eqs. (24) and (25) must vanish, 

and thus 


(-MQ 2 W +■ 


rr 


lQD rr ) (-MQ 2 


c O W 


+ e K + W COS<p„ 
0 <jxp 0 


♦ IQ 


C n 0 ) 

0 (jxp 


- (C ft + W Sin « + IQ c D ) (K + 1QD J = 0 (26) 

o r<p o o r<p <pr <pr 

Equating Imaginary and real terms of Eq. (26) to zero, we get the following 

two equations 
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M W = 


fi ‘(0 ♦ D rr ) 

w rr 


T (K-r 0 + 0 rr K ) - (B„ K „ 

+ U \ L rr w rr <p<p ' r<p q>r 


1 


+ D _K r ) - — (D W sin <p - D W cos «p )"] (27) 

<pr rip c <pr o rr o J x 


(A «) & 


4 -» 2 [»»(^^ rr * , ^ p a ) 


+ (D r J5 -DC) 
rr <p<p rip ipr 


j 


i 


+ ( K K - K K ) + (K W cos <p - K W sin ») = 0 (28) 

rr <p<p <pr r<p e Q rr o <pr o ' 


From the above two equations M and ft can be calculated. M Is the 


critical mass parameter above which the bearing Is unstable. 

Study Under Variable Viscosity Oil 

In the foregoing sections a scheme was devised for the analysis of the 
constant viscosity oil. We now consider the case of variable viscosity. The 
viscosity of most oils Increases with pressure and the following relationship 
Is assumed 

n = n Q e aP (29) 

where n Q viscosity of the oil at the Inlet condition 
a piezo-viscosity coefficient. 

The Reynolds equation under steady-state conditions for the variable 
viscosity case Is 


where 

and 

The constant 


( R 


k £ 


* i ( 
o b y 


b = 


, -p B 
1 - e p o 


n Q a <*>R 


3 

3Z / 


8 P.0 

36 


B can also be expressed as 


B = F • G (|) 


(30) 

(31) 


13 



where G 


= aE . 

It may be seen Eq. (30) Is of the same form as that of Eq. (15). In this 
case q may be referred to as the modified pressure. We can solve q Q 
following a similar approach. After knowing q Q the pressure distribution 
P q can be found using Eq. (31) for a given value of B. The elastic 
deformation 6 Is then obtained from Eq. (11). The process was repeated 
until the film shape and pressure were converged simultaneously. The 
stability study was then made using the method given In the previous section. 

RESULTS AND DISCUSSION 

The circumferential pressure distribution obtained from the present 

method of solution has been plotted In Fig. 2 for various values of F for a 

particular bearing. It may be seen that the peak pressure decreases with 

Increase In the elasticity parameter (I.e., for softer shells). In Table I 

the comparison of maximum centerline pressures In the circumferential 

direction of the present solution for a finite bearing with L/D = 1.0, 

c q = 0.85, and for values of F varying from 0 to 0.4 with those of 

Ref. (4) are shown. Table II compares the mass parameter of the present 

solution with that of Allaire (Ref. 16) for a rigid bearing. Both tables 

Indicate a fair agreement. The film thickness at the centerline In the 

circumferential direction Is plotted In Fig. 3. When F = 0.4, the film 

profile at the high-pressure zone shows a similarity to that of an 

elastohydrodynamlc point contact condition. Fig. 4 shows the steady-state 

load capacity variation with elasticity parameter for three eccentricity 

ratios (I.e., c q = 0.6, 0.8, 0.85). Although there Is little variation of 

load with F at low eccentricity ratios, the load drops very sharply with F 

at e = 0.85. This load goes below the one that Is obtained at 
o 

c q = 0.8 when F approaches a high value. The Increase In F Increases 
the minimum film thickness. This In effect reduces the true eccentricity 
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ratio, therefore pressures and load capacity drop. It may oe mentioned that a 
similar observation has been made by Conway and Lee (Ref. 6) while analyzing a 
flexible bearing using the short bearing approximation. Because of the 
difficulty In convergency. It was not possible to study the effect of load at 
values of c Q > 0.85. As there Is no effect of elastic parameter on the 
bearing performance for c q < 0.6, the results are presented for c q > 0.6. 

In Figs. 5 and 6 the stiffness and damping coefficients for the 
Isovlscous case have been plotted for e q = 0.6 and 0.8. It Is evident, 
there Is little variation of these coefficients for e 0 = 0.6. However, 
some stiffness and damping coefficients (K rr , K r<p , D rr ,) at e 0 = 0.8 
vary significantly with F. An Increase In distortion of the liner causes the 
direct-coupled radial stiffness coefficient to decrease while the direct 
damping coefficient Increases. The reverse effect was noted for the 
cross-coupled r-«p stiffness and damping coefficients. Further, all other 
stiffness and damping coefficients were relatively unaffected by the 
flexibility of the bearing liner. 

Eccentricity ratio plays an Important part In the margin of stability. 
This can be seen from Fig. 7. The stability threshold falls rapidly with F 
at higher eccentricity ratios. When one studies the stability behavior of a 
rigid bearing. It Is revealed that the stability Is not a problem at high 
eccentricity ratios. In other words, a rigid bearing Is highly stable at high 
eccentricity ratios. This can also be verified from Table II. However, for a 
flexible bearing liner the converse Is true. The whirl ratio (Fig. 8) also 
Increases very sharply with F for e Q = 0.85 and It goes well above 0.5. 

This Is a perfect example of the latitude that can be taken with the meaning 
of "half-frequency whirl" since the whirl ratio for Instability occurs at 
values generally close to, but different from 0.5. 
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The effect of L/D ratio on stability Is shown In Fig. 9 for a 
particular value of F = 0.2. It Is seen that the bearing Is highly stable 
when L/D Is small but drops drastically as L/D Increases from 0.5 to 1. 
The stability margin does not decrease much when L/D Increases from 1 to 2. 
A similar trend Is also seen for other values of F. 

In Figs. 10 and 11 the variation of stability and whirl ratio on the 
pressure-dependent viscosity Is shown. The stability Is seen to Improve with 
Increasing B for all values of e Q while the whirl ratio generally 
decreases with Increasing B. 

CONCLUSIONS 

Numerical methods are used to determine the effects of elastic 
distortions In the bearing liner on bearing performance and stability for a 
finite journal bearing. A linearization method using a first-order 
perturbation theory was used for the stability analysis. The following 
conclusions were evident: 

1. The region of stability decreases as the bearing liner Is made more 
flexible for high eccentricity ratios (l.e., c q > 0.8). For c q < 0.8, 

the flexibility of the bearing liner had little or no effect on stability. 

2. As L/D Is Increased, distortion effects are more prominent. This 
leads to a decrease In stability. 

3. An Increase In distortion of the liner causes the direct-coupled 
radial stiffness coefficient to decrease while the direct damping coefficient 
Increases. The reverse effect was noted for the cross-coupled r-<p 
stiffness and damping coefficients. Further, all other stiffness and damping 
coefficients were relatively unaffected by the flexibility of the bearing 
liner. 
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4. The hydrodynamic pressure and hence the load capacity la reduced as 
the bearing liner becomes more flexible, especially at eccentricities greater 
than 0.8. 

5. As the bearing liner becomes softer, the whirl ratio, Q, remains 
relatively constant (l.e., 0.46 < Si < 0.53) for c Q = 0.6 and 0.8. For 

c q = 0.85 the whirl ratio changed drastically (l.e., 0.30 < a < 0.58) with 
an Increase In the elasticity parameter. 

6. The stability Increases for Increasing viscosity of the oil. 
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TABLE I. - COMPARISON OF MAXIMUM 
STEADY-STATE PRESSURE OBTAINED 
BY THE PRESENT METHOD TO THAT 
OF BRIGHTON ET AL. (4) 


[L/D = 1.0, c 0 = 0.85, 
H/R = 0.3, a = 0.4] 


F 

p max 

p max a 

0 

17.25 

16.80 

0.05 

13.50 

14.10 

.10 

11.50 

11.40 

.20 

9.00 

8.70 

.40 

6.25 

6.30 


indicates the results 
by Brighton et al. (Ref. 4). 


TABLE II. - COMPARISON OF THE 
PRESENT SOLUTION WITH THAT OF 
ALLAIRE (16) FOR RIGID 
CYLINDRICAL JOURNAL 
BEARING 


[L/D 

= 1.0, 

F = 0] 

c O 

M 

M a 

0.1 

6.950 

6.8 

.3 

7.512 

6.8 

.8 

15.45 

16.0 


indicates the results 
by Allaire (Ref. 16) . 
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FIGURE 1. - SCHEMATIC DIAGRAM OF A JOURNAL BEARING. 
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FIGURE 2. - CIRCUMFERENTIAL PRESSURES AT THE CENTERLINE 
(ISOVISCOUS). 



FIGURE 3. - FILM THICKNESS AT THE CENTERLINE (ISOVISCOUS). 
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FIGURE 7. - VARIATION OF STABILITY WITH ELASTICITY 
PARAMETER (ISOVISCOUS) . 


e 0 



FIGURE 8. - VARIATION OF WHIRL RATIO WITH 
ELASTICITY PARAMETER (ISOVISCOUS). 
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L/D = 1.0, F = 0.1, 0 = 0.4, H/R = 0.3 



FIGURE 11. - EFFECT OF VISCOSITY OF OIL ON WHIRL RATIO. 
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